CHAPTER VII 

Theoretical Analysis of Modes of Vibration for Isotropic 
Rectangular Plates Having All Surfaces Free 

By H. J. McSKIMIN 

7.1. Introduction 

The comparatively recent advent of crystal controlled oscillators and of 
wave filters employing piezoelectric elements has resulted in an extensive 
study of the ways in which plates made of elastic materials such as quartz 
or rochelle salt can vibrate. Of special interest have been the resonant 
frequencies associated with these modes of motion. As will be indicated in 
subsequent paragraphs, the general solution to the problem of greatest 
interest is quite complex, and has not been forthcoming, (i.e., as applied to 
rectangular plates completely unrestrained at all boundary surfaces). For 
this reason numerous approximate solutions have been developed which 
yield useful information in spite of their limitations. Several of these 
solutions will be discussed in the following sections. The three general 
types of modes (i.e., the extensional, shear, and flexural) will be analyzed in 
some detail. Also, as a preliminary step the formulation of the general 
problem along classical lines will be developed. 

For the most part, the solutions obtained here are limited to those for an 
isotropic body. However, such solutions provide considerable guidance 
for the modes of motion existing in an aeolotropic body such as quartz. 

7.2. Method of Analysis 

In order to set up the desired mathematical statement of our problem it 
will be necessary to consider first of all two very fundamental relationships. 
The first of these is the well known law of Newton which states that a 
force/ acting on a mass m produces an acceleration a in accordance with the 
formula 

/ = m-a 

The second relationship which we shall need is Hooke's law relating the 
strains in a body to the stresses. If forces are applied to the ends of a long 
slender rod made of an elastic material such as steel (Fig. 7.1) a certain 
amount of stretching takes place. If the forces are not too great, a linear 
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relationship between the applied stress and ensuing strain is found to exist. 
Expressed as an equation 

v 

— = E in which X x is the force per unit area, 
x x 

x x is the strain per unit length, and £ is a constant known as Young's 
Modulus. (Refer to Section 7.7 for further definition of terms). 

In an analogous manner, shearing stresses applied to an elastic solid as 
shown by Fig. 7.2 produce a shearing strain such that 

— = A, the shear modulus. 
x v 

In general there will be contributions to a particular strain from any of 
the stresses which may happen to exist. For example, when an isotropic 



Fig. 7.1 — Bar under tensional stress 




Fig. 7.2 — Bar under shearing stress 

bar is stretched, there will be a contraction along the width which has been 
produced by a stress along the length. A statement of these relationships 
(known as Hooke's Law) is given by the equations of Section 7.8. 

It is now of interest to consider the conditions of equilibrium for a very 
small cube cut out of the elastic medium which in general is stressed and in 
motion. Reference to Fig. 7.3 will help to visualize the stresses which may 
exist on the faces of this cube. Since these stresses vary continuously 
within the medium, a summation of the forces acting on the cube along each 
of the major axes can be made with the use of differential calculus. From 
Newton's Law previously cited, it is apparent that any unbalance of these 
forces will result in an acceleration inversely proportional to the mass of 
our small cube. Three equations may then be derived, one for each major 
direction. 1 If only simple harmonic motion is considered (i.e. all displace- 

1 Refer to "Theory of Elasticity" by S. Timoshenko or to any standard text on elasticity. 
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ments are proportional to sin col where to = 2t times frequency) the following 
simplified equations result. 




Fig. 7.3 — Stresses acting on small cube 



dX x _|_ dX v . dX t _ 
dx dy dz 

dy dx dz 
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dZ t dX. 3 Y. 
dz dx dy 



= —pio w 



(7.1) 



Since stresses are related to strains in a very definite manner, the above 
equations may be converted into a more useful form involving only displace- 
ments. For isolrQpic media, the following results. 
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dx 



— pco u 



A^v -\-B^ = - P uv 
dy 

AV 2 w + B — = —paw 
dz 



(7.2) 



In this grouping, 



6a; 2 T fly 2 dz 2 

cte dy dz 

and A and 5 are given in terms of the fundamental elastic constants X and 
p with A = ju, B = X + p. 

An even more elegant statement of the equilibrium conditions attributable 
to Love 2 follows immediately from equations 7.2, since by differentiating 
each one in turn with respect to x, y, and z respectively, and then adding 
results, one obtains the wave equation 



(V 2 + h*)e - o 



(7.3) 



where 



h 2 = 



pd) 



pOi 



A + B X + 2m 



Whatever our solution may be, then, it must satisfy equation (7.3). If 
such an expression for e is found, the displacements formed in the following 
way will satisfy equations 7.2 as can be shown by direct substitution. 



U h 2 dx 



_1 de 

V " W dy 



1 de 

w= -ltTz 



(7.4) 



In addition to equations 7.2, another set of requirements will be necessary 
when any particular problem is considered. They are known as the bound- 
ary conditions, and in general are easily deduced from a knowledge of how 
the plate or bar is held. 

For a rectangular plate free on all surfaces, the boundary condition is 
simply that all surface tractions vanish. This requires certain stresses to 
become zero at the boundary as can be seen from the following expressions 
for the x, y, and z components of traction in terms of unit stresses. 

1 A E. Love, "A Treatise on the Mathematical Theory of Elasticity." 
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X = X z t + X h m + X t n } 
Y = Y y m + Y t n + XJ 
Z = Z t n -f X t e + Y t m 



> = for free surfaces 



(7.5) 



((, m, and » are direction cosines of the normal to the surface at the point in 
question). 

The general problem is now seen to be one of finding solutions for the 
displacements u, v, and w such that both the equilibrium and boundary 
conditions are satisfied. In the following section several interesting solu- 
tions will be considered for rectangular plates having all surfaces free, this 
being the case of greatest interest in so far as this paper will be concerned. 

7.3. EXTENSIONAL VIBRATIONS 

One of the most useful modes of vibration of practical interest is the 
extensional, in which particle motion takes place in essentially one direction 
so as to alternately stretch and compress the elastic medium. Piezoelectric 




x = l 



Fig. 7.4 — Longitudinal bar 



plates vibrating in this manner, and of the shapes shown in figures 7.4 and 
7.5 have been used extensively in wave filter and oscillator circuits. The 
approximate resonant frequencies corresponding to this type of motion are 
easily obtained by a consideration of equations 7.1 and 7.2. For the 
longitudinal bar of Fig. 7.4 the only stress that need be considered is the X x 
extensional, all other stresses being so small that they can be neglected. 
The equilibrium equation then becomes 



or, since 



ax x _ 

dx 


2 

= — pu U 


du _ 

dx = 


±X 


dx 2 


P 2 „. 
--uu 



(7.6) 



(7.7) 
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It is easily seen that u = cos kx is a solution to this equation if k = w Kf % • 

If now the boundary condition that the stress X x must become zero at the 
ends of the bar (i.e., x = 0, x = t - refer to Eq. (7.5)), is fulfilled, the 

solution will be complete. At % = 0, X x = E — will always equal zero. 

Furthermore, \i k — ^ or any whole number multiple of - the extensional 

stress will likewise reduce to zero at x = I. The desired solution will then 
be as follows, / being the resonant frequencies. 



u = COS O) 



m = 1, 2, 3, etc. 



(7.8) 



•HtK- 




-*-x 



Fig. 7.5— Thin plate 

The plate of Fig. 7.5 will now be considered. Here it can no longer be 
assumed that the X x stress is the only one of importance. Instead, the 
displacements v and w will be considered zero and the displacement u a 
function of x only. This means that the shear stresses X u , X„ , Y t vanish, 
so that the equilibrium equations 7.2 reduce to 



d 2 u 



^2 

o u 



or 



.D» . n U» 2 



d 2 M — pw U 



(7.9) 



(7.10) 



dx 2 A + B 

This is seen to be of the same form as equation (7.7) previo usly dis cussed, 
and will again have the solution u = cos kx with k = u A/ 



The 
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boundary condition on the X x stress will be met if k = — so that the 

t 

following solutions result. 3 



/p(l +<r)(l - 2a) J p~~ 

U = C0SW V E(l-<r) * = C ° S " V JT+1S* 

n f _ ™* t /E (1 - *) mw /X + 2 M (7.11) 

» - 2x/ - T y -- (1 + (r)(1 _ 2(r) = T j/ -y- 

m = lj 2, 3, etc. 

It is seen that this formula for resonant frequencies is the same as given by 

equations 7.8, with £ replaced by ^ — r , so that the frequency 

(1 — 2<r* — a) 

constant/-/ will be somewhat higher than f-l for a long slender bar. 

It is recognized that the solutions derived above hold true only for the 
limiting cases of a long slender bar, and a very thin plate respectively. It is 
therefore of interest to trace the resonant frequencies corresponding to these 
extensional modes of vibration as departure is made from the limiting cases 
mentioned above. 

An experimental plot of the resonant frequencies of a thin plate of length I 
and width w reveals that the frequency of the longitudinal mode first 
discussed is gradually lowered as the width of the plate is increased. There 
is also another frequency corresponding to an extensional vibration along 
the width which for a very narrow plate corresponds to the second type of 
extensional mode considered in the foregoing paragraphs, except that the 
frequency constant will be slightly different because coplanar stresses are 
involved. 

As seen from Fig. 7.6, the resonant frequency curves do not cross, but 
exhibit coupling effects. This is understandable from the fact that motion 
in one direction is mechanically coupled to motion in the other as indicated 
by Poisson's ratio a. 

In order to derive expressions for the u and v displacements associated 
with the extensional mode along the length, taking into account the above 
coupling effect, the following analysis proves interesting. 

Consider the infinite isotropic strip of width b as shown by figure 7.7. 
As will be demonstrated presently, solutions can be found such that the 
equilibrium equations and the boundary conditions are precisely satisfied. 
Furthermore it will be found possible to cut a section out of this strip in 

3 If the length and width of the plate are very large in comparison to the thickness, 
the boundary conditions for the Y v and Z, stresses may be neglected without causing 
appreciable error. The quantity A + B has been evaluated in terms of E and a for 
purposes of comparison. 
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such a way that the boundary conditions for the cut edges are very nearly 
satisfied. The plate formed in this way may then be considered as vibrating 
at the required frequency /, which will then be the resonant frequency 
desired. 
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Fig. 7.6— Extensional modes with mechanical coupling 
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Fig. 7.7 — Infinite strip 



Let displacements be arbitrarily chosen in the following way: 
u = U cos kx cos („ 
v = V sin kx sin /„ . 



(7.12) 
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As shown in Section 7.9, two solutions of this type will satisfy the equi- 
librium equations precisely. One corresponds to e = in the wave equation 
7.3, while for the other e j* 0. Superposition of these two solutions and 
proper evaluation of parameters make it possible to satisfy the boundary 

conditions at the edge of the strip; namely, that at y = ± -, F„ = and 

X y = 0. (Refer to equations 7.5). The following transcendental equation 
is obtained 



, t b 
COth 2 -2t x ttk\\-v) 



in which 



cot /, \ <*-«<« + •« 



«-m'-* 



(7.13) 



a = f-j 



2 



(7.14) 



This equation may be solved graphically to yield values of frequency 
corresponding to given values of k. For our discussion of the length ex- 
tensional mode of vibration, the first root only will be considered. 

Fig. 7.8 shows a plot of 6-b against b-k assuming that Poisson's ratio is 

.33. 4 If k = 1, and b = 1, for example, 6 = a/-? w = 1.62. 

The equations for the displacements when determined as explained in 
Section 7.9 become: 

u = Vi [cosh -344 y + .402 cos 1.278 y] cos x 

(7.15) 
v = Ut [.344 sinh .344 y + .315 sin 1.278 y] sin x 

All three stresses X X) Y v , and X v may be calculated from the above 
equations. If the length of our plate is made equal to mir, where m is an 
integer, the extensional stress X x will equal zero regardless of y at the 
boundaries x = and x = I since X z a sin x = when x = mir. Also it 
can be shown by calculation that X y is so small in comparison to the exten- 
sional stresses as to be entirely negligible; hence our solution is complete. 

If I = x, the plate will be vibrating in its fundamental longitudinal mode. 
The distortion which results is shown by Fig. 7.9. It is seen that most of 

4 Plotted in this way, the same curve results regardless of the value of b chosen for the 
purpose of solving Eq. 7.13. 
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the motion is along the x axis, though there is a certain amount of lateral 
contraction as the plate elongates. 

The second harmonic will have the same resonant frequency if t = 2ir, 
the third if t = 3ir, etc. 
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Fig. 7.8 — 0-b versus k-b for plate longitudinal modes f a = .33, & = 
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Fig. 7.9 — Distortion of plate vibrating in first longitudinal mode 

In addition to harmonic modes along the length just considered there 
will be those for which the motion breaks up along the width. In general, 
the distortion of the plate may be quite complex with simultaneous variations 
along both dimensions. Similarly, for plates such as shown in Fig. 7.5 
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there will be many extensional modes which have resonant frequencies 
somewhat above those given by Eq. 7.11. Analysis of the motion shows 
that for these modes the displacement along the thickness varies periodically 
(or "breaks up") along the major dimensions of the plate. There again the 
distortion pattern of the plate may become very complex. 

7.4. Shear Vibrations 

The second class of vibrations which will now be considered is the shear. 
This type of mode is of special importance because of the fact that piezo- 
electric plates vibrating in shear are widely used for frequency control of 
oscillators. For example, the AT quartz plate which is so much in demand 
utilizes a fundamental thickness shear mode in which particle motion is 
principally at right angles to the thickness. The distortion of the plate will 
be similar to that shown in Fig. 7.2. 

A simple, yet very useful formula for the resonant frequencies associated 
with the above type of displacement has been derived on the assumption 




Fig. 7.10 — Orientation of thin plate 

that the length and width of the plate are very large in comparison to the 
thickness. For the xy shear mode, the displacement u is assumed to be 
u = U cos ky, all other displacements being equal to zero. The only stress 
that need be considered then, is the X y shear which is proportional to sin ky. 
Boundary conditions on this stress at the major surfaces of the plate are 
easily satisfied by choosing k such that X y = at y = and y = t. (Refer 

to Fig. 7.10.) This will be the case if k = —-, where m is any integer, and 

t is the thickness of the plate. By using the simplified equilibrium equation 
as reduced from equations 7.1, a formula for the resonant frequencies is 
obtained in much the same manner as for extensional thickness modes. 

u = 2-rf = ™ J^- tn = 1, 2, 3, etc. (7.16) 

In this formula the shear modulus A appears instead of Young's modulus 
as in the case of longitudinal modes. Harmonic modes are given by values 
of m greater than unity. 
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In addition to the resonant frequencies predicted by the foregoing analysis, 
there will be others corresponding to shear vibrations in which the principal 
shear stress varies periodically along the length and width of the plate. 
A formula which yields the approximate frequencies for these modes is 
developed in Section 7.9. It is shown that if the length and width are 
large in comparison to the thickness, the following expression may be used: 

w = 2t/= *y-y ciijf+ -y- + -y (7.n) 

In this formula which has been derived for xy shears the c constants are 
the standard elastic constants for aeolotropic media. For isotropic plates 
such as have been considered up to this point 

„-*$£•*.. x + * 

1 — 2a* — a 

and 

Cf& = cee = A, the shear modulus (7.18) 

Various combinations of the integers m, n, and p may be chosen, with the 
restriction that neither m nor n can equal zero. It is seen that if t and w 
are very large the formula reduces to that of Eq. 7.16 which was derived on 
precisely that basis. Also, it is seen that the more complex modes all lie 
somewhat above the fundamental shear obtained by setting m = n = 1 
and p = 0. 

Plate shear modes are also of considerable interest, particularly the one 
of lowest order. For a plate having a large ratio of length to width a formula 
similar to that given by equation 7.17 (but for two dimensions only) may 
be developed. If the plate is nearly square, however, this formula no longer 
yields sufficiently accurate values for the resonant frequencies. Coupling 
to other modes of motion 5 complicates the problem so much that only 
experimental results have been of much practical consequence. Fig. 7.11 
shows in an exaggerated way the distortion of a nearly square plate vibrating 
in the first shear mode. 

7.5. Flexural Vibrations 
7.51. Plate Flexures 

One of the most studied types of vibrations has been the flexural. Perhaps 
this is true because it is the most apparent and comes within the realm of 
experience of nearly everyone. The phenomena of vibrating reeds, xylo- 
phone bars, door bell chimes, tuning forks, etc. are quite well known. 

8 It is found experimentally that odd order shears are strongly coupled to even order 
flexures; similarly, even order shears and odd order flexures are coupled. 
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Beam theory has been used quite extensively to derive the equations 
which yield the resonant frequencies and displacements for bars vibrating 
in flexure. To obtain reasonably accurate results for ratios of width to 
length approaching unity, however, the effects of lateral contraction, rotary 
inertia, and shearing forces must be considered. This leads to a rather 
complicated solution which is much more accurate than that derived by the 
use of simple beam theory only, though it is still approximate in nature. 

For two dimensional plates free on all edges a method of analysis may be 
used which is similar to that described under extensional modes. While 
it is somewhat involved it yields direct expressions for the two displacements 
u and v, so that all stresses may be calculated, and the extent to which 
boundary conditions are satisfied determined. 




Fig. 7.11 — Distortion of plate in first shear mode 

Solutions for u and v are assumed to be of the form 

u = U sin (y ccs kx 

(7.19) 
v = V cos fy sin kx 

For the infinite strip previously considered a transcendental equation is 
obtained which is the same as equation 7.13 with the exception that the 
left-hand expression is inverted. 



tan/l 2 _ -2f^,k\l - a) 

t ,b (d- k 2 )(Cl + ok 2 ) 
tan & -z 



(7.20) 



(Refer to Eq. 7.14 also.) 

8 This is an extension of DoerfTler's analysis used to obtain harmonic flexure frequenceis 
for plates — "Bent and Transverse Oscillations of Piezo-Electxically Excited Quartz 
Plates"— Zeitschrift Fur Physik, v. 63, July 7, 1930, p. 30. Also refer to "The Distribu- 
tion of Stress and Strain for Rectangular Isotropic Plates Vibrating in Normal Modes of 
Flexures"— New York Univ. Thesis by Author, June 1940. 
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The lowest order solution to this equation is found to correspond to 
flexure vibrations in the infinite strip. A calculation of stresses, however, 
reveals that boundary conditions cannot be satisfied properly even for the 
case of a long narrow plate. It can be shown, however, that another 
solution may be derived for the same value of frequency by letting k become 
imaginary. This simply means that the u and v displacements become 
hyperbolic functions of x instead of sinusoidal. The two complete solutions 
for the infinite strip may then be superimposed and parameters adjusted so 
that for definite values of length corresponding to fundamental and harmonic 
modes the proper stresses reduce essentially to zero on the ends of the plate. 
For plates having a ratio of width to length less than .5, this method gives 
very accurate expressions for displacements and stresses. If only the 
resonant frequency is required, ratios up to unity and beyond (for the 
fundamental mode) may be considered. 

An example has been worked out to provide a complete picture of the 
displacements for a bar of width = 1, k = 1 and a = .33. Use of equation 

7.20 yields the quantity 6- = — = .166 from which the resonant frequency 

may be obtained. Using this value of 2 , one finds that k 2 = —.800 also 
satisfies equation 7.20. By making the total length of the bar equal to 
4.50 the X x extensional stress and the X y shear stress may be made essen- 
tially zero on the ends of the plate regardless of y. 
The following expressions for u and v are obtained: 

u = (sinh .9132 y - 1.02 sinh .9718y) sin x 

-.160 (sin .9828y - .9568 sin .9250y) sinh .8944a; 

(7.21) 
v = (-1.094 cosh .9132y + .9915 cosh .9718) cos x 

-.160 (.9095 cos .9828y - .990 cos .9250y) cosh .8944a; 

Fig. 7.12 shows the distortion of the plate as calculated from the above 
expressions. It is seen that there will be two points at which there is no 
motion in either the x or y directions. These nodal points can be used in 
holding the plate, since it may be clamped firmly there without altering 
the displacements or resonant frequency. For the example shown, these 
nodes are positioned a distance of .21 1^ from the ends of the plate as com- 
pared to .224f for a long thin bar. 

7 A graphical solution to determine I is most convenient in which parameters are 

adjusted so that X x = at x = ±~ and y = ±^; X u = at x = ±- and y = 0. These 

w 
stresses will remain essentially zero for all values of y if the ratio of j is not too great. 



MODES OF VIBRATION 



165 



Figures 7.13 and 7.14 show the distribution of the principle stresses as a 
function of position along the length. It is seen that for the particular 




DISPLACEMENT SCALE 

Fig. 7.12 — Distortion of bar vibrating in first free-free flexure mode 
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Fig. 7.13— Distribution of longitudinal stress for free-free bar vibrating in first 

flexure mode 



T* 



o 

0.001 
0.002 



>■ 0.003 
0.004 
0.005 

a006 o^ o^ 06 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

X 

Fig. 7.14 — Distribution of shear stress for free-free bar vibrating in first flexure mode 

example cited, the maximum shear stress is only about one-tenth the 
maximum X x extensional stress. Both of these stresses reduce to zero at 
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the ends of the plate as they should in order to satisfy the boundary condi- 



w 
tions. As the ratio of — is increased the shear stress becomes of greater 

importance. 

7.52. Thickness Flexures 

The final analysis to be considered in this paper is for thickness flexures 
along the width or length of a thin plate. These modes are of particular 
interest in connection with the dimensioning of quartz plates for which it is 
desirable to utilize the fundamental thickness shear mode. (AT plate, for 
example.) It is found experimentally that even ordered thickness flexures 
are coupled to this shear to such a degree that at certain ratios of dimensions 
the operation of the plate as an oscillator or filter component is impaired. 

The two-dimensional solution derived in the preceding paragraphs can 
be used to predict certain harmonic thickness flexures; however, in order to 
obtain a complete picture it is necessary to extend the theory to three 
dimensions. This has been done by the author with the following transcen- 
dental equation as a result (refer to Section 7.93). 



tan t\ - 

It 



-llxhAo? 



tan /« 



b UBifi + a 2 ) + Ami ~ a] 



(7.22) 



Solutions to this equation are exact in nature for a plate of thickness b 
and of infinite extent in both the x and z directions. The quantity a 2 is 
equal to the sum of the squares of k and m which appear in the expressions 
for displacements as follows: 



(7.23) 



u = V fi(y) sin kx cos mz 


v = V fiiy) cos kx cos mz 


w = W fi(y) cos kx sin mz J 


Also in equation (7.22) 


f 2 Q 1 A 2 




(I = 6 2 -a 2 


> 


A 





(7.24) 



The lowest order solution to equation (7.22) with a 2 positive again cor- 
responds to flexure vibrations, as in the two dimensional case. Fig. 7.15 
shows a plot of 8- b against a-b calculated for a = .3. 
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For reasonably high order flexures it may be reasoned that the true dis- 
placements will be very nearly the same as those for the doubly infinite 
plate as derived by the above method since the correction necessary to 
fulfill the boundary conditions will only apply very close to the edges of the 
plate. It will then be sufficient to choose values for k and m such that 

k = p - and m = — where p and q are integers. The values of a 2 obtained 

r I w 

in this way determine the corresponding resonant frequencies. 



























































































































T 


























































n 






































































/ 


/ 


























/ 


/ 

















































































<xb-bVk 2 +m 2 
Fig. 7.15 — 0-b versus a-b for thickness flexures 

If it is desired to solve for the ordinal xy flexures, for example, m should 
be set equal to zero. The displacements in this case will be independent 
of the z dimension. When q is assigned values other than zero however, 
the resulting modes may be considered as xy flexures which vary or break 
up along the third major dimension. If q is small the resonant frequencies 
will lie only slightly higher than that of the corresponding ordinal flexure 
for which q = 0. 

Fig. 7.16 shows a few of the resonant frequencies as calculated for values 
of shear modulus and density corresponding to AT quartz. The effects of 
coupling to the fundamental thickness shear are shown by dotted lines for 
the 14th xy flexure. As might be expected there is similar coupling between 
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the 14th flexure which breaks up once along the z dimension and the shear 
which breaks up once along z — etc. A few of these flexures which break up 
along 2 are shown for the 16th ordinal flexure. 
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Fig. 7.16 — XY thickness flexure modes for square plate 



7.6. Summary 

Three main classes or families of vibrational modes are found to exist in 
rectangular elastic plates free on all surfaces; namely, the extensional, the 
shear, and the flexural. In general, the associated displacements are 
functions of all three dimensions and may vary in such a manner as to make 
the distortion of such plates quite complex. 

For certain limiting cases, approximate solutions for the resonant fre- 
quencies and displacements (from which strains and stresses may be cal- 
culated) can be derived. Though there are a number of methods that can 
be used for specific problems, it has been found very convenient to utilize 
the classical formulation. For this reason the basis of this method has been 
discussed briefly. In essence it requires that displacements and stresses 
occurring within the elastic solid satisfy conditions of equilibrium as de- 
rived from Newton's Law. At the boundaries, certain other relations must 
be satisfied in order that conditions of clamping might be fulfilled. For 
plates entirely unrestrained the latter requires that all forces (tractions) 
acting through the free surfaces must vanish. 

For thin rectangular plates (such as quartz crystal oscillator plates) the 
modes of greatest practical consequence are plate modes, for which all 
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stresses are essentially coplauar and independent of the thickness, and 
thickness modes, for which all dimensions must be considered except in 
limiting cases. 

Because of their great utility, simplified formulae have been derived for 
the resonant frequencies associated with long, narrow bars vibrating longi- 
tudinally, thin plates with extensional motion along the thickness dimension, 
and thin plates vibrating with shearing motion at right angles to the 
thickness. 

Exact solutions for the infinite strip have been derived, and used in 
obtaining the displacements and resonant frequencies for flexural and 
longitudinal modes. Such solutions take account of the fact that the width 
of the plate may become appreciable. While limiting cases of plate shear 

may be analyzed, solutions for ratios of — approaching unity have not 

V 

proved very satisfactory. This is attributable to the fact that coupling to 
flexural modes is severe. 

Thickness flexural modes which exhibit displacement variations along 
both length and width dimensions of the plate have been analyzed by 
extending the "infinite strip" theory to three dimensions. Solutions 
obtained are fairly accurate if the harmonic order of the flexure is sufficiently 
great. 



7.7. Nomenclature 



p = density 
E = Young's modulus 
a = Poisson's ratio 

A = Shear modulus = 

E 



E 



= M 



B = 



2(1 + cr)(l - 2<r) 
E 



2(1 + <r) 

= X + n for 3 dimensions 



for plane stress 
2(1 - a) 

co = angular velocity = 27r/ 



,2 pw_ 
A 



u, v, w = displacements in x, y and s directions 

du . dv . dw 
dx dy dz 
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V = Laplacian =_ + _ + - 

fV> AJ Z \ 

1 ' ' > unit strain components 

*r, *., yzj 

* ' "' * \ unit stresses 
X v , X, , Y t J 

7.8. Stress-Strain Equations for Isotropic Media 
Xx = ^ (X* - <rF y — <rZ,) 

yv = 2 ( Y v — ffX * ~ aZ ') 

z t = •= (Z, — o-X* — oT„) 

XL 

shear strain = - X shear stress 
A 



x x 


->(* + 'S 


Yy 


- a (-*+'S) 


z t 


-«(*+'S 


Xy 


-<g+i) 



*-'(* + =) 

F, = A (- + — ^ 
\dz dy/ 



For plane stress in xy plane 



*» « i (X x - oY v ) 

y u = jj (^v - ^x) 
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x z 


E 
l-o 2 


(x x 


+ <ry v ) 


Yy 
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1 _2 


(y. 


+ <TX X ) 
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7.9. Mathematical Derivations 
7.91. Longitudinal Vibrations in Two-Dimensional Plates 
As explained in the text, solutions for the infinite strip of Fig. 7.7 are 



first derived. Let 



u = U cos kx cos l v 
v = V sin kx sin £„ 



(7.12) 



_ , . du , dv 

where U and 7 are constant. From these expressions e = — + — can 

ox ay 

be obtained and substituted into the equilibrium equations 7.2. Two 



expressions as follows result after dividing 
cos kx cos l v . 

Bk 



A{k* +e)-i£ (-kU - IV) = pu 
A(k 2 +/) + 5f (-JH7 + /7) = pco 



Subtracting the second from the first of these equations, it is seen that 



# + f)<- 



kU + 



rr) =o 



Either or both of these factors equal to ze o will satisfy 7.26, so that two 
V 



I back into equations (7.25), 
conditions on co 2 are found. The two solutioi is will be 



values of — are obtained. By substitutin 



U\ k 



t/2 €z 



By superimposing the two solutions the u ai id v displacements now become 
u = [Ui cos l\ y + C/2 co 1 4 y] cos &.? 

r u k 1 ( 7 - 28 ) 

a = — - Ui sin Ay + y #2 sin / 2 y sin ftg 



hrough by the common term 



(7.25) 



(7.26) 



+ t\) = 



pu 



(7.27) 



. 



• 
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Using the relationships of Section 7.8, one may now calculate all stresses. 
The argument k has purposely been kept the same for both of the super- 
imposed solutions in order that boundary conditions at y = ±- might be 

satisfied regardless of x. 

For Y v to equal zero at the edges of the strip 



dv , du 
dy dx 

This gives rise to the equation: 



^6 = (7.29) 



J7i(/f + ak 2 ) cos it \ - U2[k\l - a)] cos & ^ = (7.30) 

— z 

Similarly, if X y = at y = ±- 



du , dv 
dy dx 



,h = (7.31) 

V 2 



Another relation is obtained from (7.31): 

-2/i/ a Ui sin l x - + U 2 (k 2 - t\) sin ( 2 h - = (7.32) 

— z 

The two equations (7.30) and (7.32) will be satisfied if the determinant 
of the coefficients of Ui and Z7 2 vanishes. The following transcendental 
equation will then be obtained; values of A 2 and 4 2 being those required 
by Eq. 7.27. 

,, b 

COt *2_ -2Mk\l - a) 

cot A ^ 

By using either of equations (7.30) and (7.32), one may derive the relation 
between Ui and TJ\ provided a solution to (7.13) is found. 

—2l\ti sin t\- 

Ui = Ui ?- (7.33) 

{t\ — k 2 ) sin t% - 

To solve equation (7.13) assume a value for k 2 and plot graphically the 

2 

right and left hand expressions as functions of 2 = — . Roots are indicated 
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by the crossover points. Values of d 2 corresponding to different values of 
k 2 may also be found in this way and a curve plotted for 2 versus k 2 , (or for 
6-b versus k-b). 

7.92. Thickness Shear Vibrations 

To obtain a formula for the approximate resonant frequencies o*f thickness 

W L 

shear for a plate having large ratios of — and — , one may consider the 

following displacements: 

u = V sin kx cos l y cos rz 

v = (7.34) 

10 = 

If there are no cross couplings between shear stresses or between shear 
and extensional stresses one may write: 8 

du . dv . dw ,,. , . 

X x = en ■=- + cu t- + en — = en kU cos &* cos ly cos rz 
3a; dy az 

X« = cee ( — H ) = — Gee IU sin ^a; sin /y cos rz (7.35) 

X t = c B b(^ + -^') = — gb6 rU sin ** cos /y sin rz 
\dz dx/ 

Substituting into the first of the equilibrium equations (7.1) and dividing 
through by common factors 

en k 2 + c e6 P + g 6B r 2 = pw 2 (7.36) 

The other two of equations (7.1) may be neglected if k and r are quite 

small so that it will only be necessary to consider equation (7.36) which 

can be solved for w 2 . It will be noticed that the X„ shear stress will pre- 

, . . . . , t . ,. » nir . vn/K 

dominate under these restrictions on k and r. Letting k = — , I = — , 

and r = — («, m and p are integers) in order to satisfy the boundary condi- 
W 

tions for this stress and also for X z , one obtains the following formula. 
(This choice of k, I, and r is also required if the shear stress is to vary in 
essentially the same manner as is experimentally observed.) 

u = 2irf = Tr^-A/-jy-)r-^+y^ (7.17) 

in which L and W must be much larger than T. 

8 Refer to equation (7.18) for values of c constants for isotropic case. 
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The boundary condition for the extensional stresses will not be met; 
however, they will be quite small in comparison to X v if k is small, and may 
be neglected. 

7.93. Thickness Flexures 

Consider a three dimensional plate having a thickness b lying along the y 
direction. The following displacements are found to be of a form that can 
be made to satisfy the equilibrium equations 7.2. 

u = U sin kx sin ly cos mz ] 

v = V cos kx cos ly cos mz > (7.37) 

w = W cos kx sin ly sin mz 

Performing the operations indicated and substituting into the equilibrium 
equations give the following result: 

A(k 2 + e + *»*) + f? (kU -IV + mW) = pco 2 



U 
A(k 2 + I 2 + m) - y (kU - IV + mW) = pco 2 

A{k 2 + l 2 + m 2 ) + ?£ (kU -IV + mW) = pa, 2 
W 



(7.38) 



Subtract the second and third equations of (7.38) from the first: 

. Bk.Bt ft . Bk Bm 

th en tf + T = ° and T/-^ = ° 



or 



F I , W m 

_ = -_ an d _=- 



(7.39) 



Putting these values back into 7.38, it is seen that the following relation- 
ship must be satisfied. 



(A + B) (k 2 + l* + m 2 ) = pco 2 



(7.40) 



VI. W 

Letting — • = — - as in (7.39), another value for — may be obtained. 
U k U 

W 
The first and second equations of (7.38) will be satisfied for any ratio of — , 

so the 3rd equation is used. 



A(k* + f + »*) -f B[V + f + mkjj) m p 



u) = 



(7.41) 
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Solving for — , using (7.41) and first of equations (7.38) 

£ « and £-=*^ (7.42) 

U k U mk 

The first ratio is the same as (7.39). For the second solution to the equilib- 
rium equations, then the following relationships exist. 

I - -< and I = =^ (7.43) 

U k U mk 

When the above are substituted back into (7.38), it is found that 9 

A(k 2 + l*+ m 2 ) = pa. 2 (7.44) 



W m . • 

In a similar way, using — = — the following are obtained: 



U k 
W m V k 2 + m 



(7.45) 



U k ' U H 

with A(k 2 + P 4- w 2 ) = pw 2 as before. This is the second solution for 
e = 0. 

The three different solutions may now be combined or superimposed to 
give 
U = [Ui sin Ay + U z sin Ay + U 3 sin Ay] sin kx cos mz 

w = I — Ui j cos Ay — Uz t- cos Ay 

H — /T cos ^»y cos ^^ cos mz n 46") 
» - [l7i j sin Ay - U* ik *+ k ® sin hy 



sin Ay I 



+ Uz j sin Ay | cos kx sin wz 



In the above equations t\ = C» because of the double requirement 
of 7.44. 10 

It is now possible to calculate the stresses existing at any point. It is 
desired to choose U x , U% , and U 3 in such a manner that the boundary 
conditions at the two major surfaces of the plate are satisfied. By using 
the relations given in Section 7.8, the extensional stress F„ , and the two 
shear stresses A'„ and Y z are calculated with the use of 7.46. They are then 

9 It should also be noticed that e = f or this solution. 
io in = k 2 = kt = k 

VI L = Ml 2 = MI 3 = ttt 
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set to zero at the faces of the plate; i.e. at y = ±- . Three equations, after 
simplifying, result. 

For Y v = at y = ±- 

Ui \oB(k 2 + t\ + m 2 ) sin ( x | + A(\ sin d |1 



+ U 2 \a(\ sin A 51 - Uz \a{1? + m 2 ) sin ft U = 



(7.47) 



For X y = at y = ± - 



|2ft cos ft ? + r_/ 2 Uft cos ft 51 



(7.48) 



For F, = at y = ±- 

£1 



Ux \2hm cos a y + 1/ 2 r^ _ ^ (r 2 2 + * 2 )^ cos 6 y 



+ tfj 



[Yftro - j s (k 2 + m 2 )) cos ft ?1 - 



(7.49) 



In order for these three equations to be satisfied simultaneously a neces- 
sary condition is that the third order determinant formed by the coefficients 
of the U's vanish. That is, 



(<rB(k 2 + (1 + m) + All) sin ft £ Al\ sin ft 



-[ 



lL/» • / b 1 
1 Ai 2 sm/2- 





A(k* + *»*) sin ft 



[v, cos 4 y 



[2^ cos a y 

[('•-(^o)-a 

[2ft m cos ft |1 [Yft m - 5 (tf + * 2 )) cos ft fl 

[(/.«- ^ (^ + «')) cos r,y 



= (7.50) 
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By dividing row 1 by cos A - , row 2 by cos A ■= ; row 3 by cos A ^ and 

by subtracting the elements of row 3 from those of row 2, considerable 
simplification results. 

The full expansion gives the following equation, after further simplifying 
operations are performed. 

\(aB(k 2 + A + n?) + A(\) tan ft|l-[4 - *' - »1 

r (7 ' 51) 

+ 2AAU0fe a + w 2 ) tan a|J - 

It should be noticed that A has dropped out entirely. Actually A = A 
as previously explained. Also the expression 



[«. + £«! + *>] 



must not be zero, for in simplifying equation (7.51) it was used as a divisor. 
Equation (7.51) may be rewritten to give 



tanA 2 -2AM(* S + * 8 ) 



(7.52) 



, b [<TB(k 2 + t\ + ™ 2 ) + ilfflW - * 2 ~ *» 1 
tan A j 

In the above 

(i4 + 5) (*> + d + m«-) = /x, 2 j 

i4(jfe 2 + /5+W 2 ) = P" 2 J 
2 

By letting 2 = — and k- + m- = a 2 equations (7.52) and (7.53) above 
become 

t b 
tan fi - nl » A 2 

2 = —KihAoL , 7 22) 

~Tft = [aB{(l + a) + iUjJ'Wi - a 2 ] 

tan A 2 

/2 fl 2 J ^ 2 

with A = e ._— - a 

^ = / 3 2 = 2 - a 2 

Equation (7.22) represents the general solution for normal thickness 
vibrations in an isotropic plate of finite thickness extending to infinity in 
both major directions. The analogy for plates of finite dimensions is 
considered in the text. 
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